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FREQUENCY SEPARATION METHOD FOR RELAXATION PROBLEMS 
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The modern era in spectral line broadening began with the understanding that the 
slow(quasistatic) ion and fast(impact) electron perturbers could be treated separately. The problem 
remained ol unifying these two theoretical limits. A scheme for this unification is presented here 
that has at its foundation a fundamental observation that is supported by analytical theory and is 
further demonstrated by computer simulation. The fundamental observation is that the ions and 
electrons can be separated most of the time, and that a frequency separation within each perturber 
subsystem can be used for unification. That is, the rigorous inclusion of slow, but not necessarily 
static ions together with the correct impact ion perturbations, will produce valid ionic line shapes. 
We show that a frequency separation may be effected to exactly include the fast modulation limit in 
a variety of modern methods that can deal with the intermediate regime between the fast and slow 
frequency limits of the perturbation. 
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A large variety of systems with relaxation behavior 
jl],D can be described in terms of an active subsystem 
in interaction with an external "heat bath" HQ]. The 
interaction with the heat bath is usually taken to be a 
stochastic, time-dependent potential V(t) that does not 
depend on the subsystem of interest. This interaction 
causes a decoherence, or memory loss, described by a 
set of time evolution operators that have, after a char- 
acteristic time, essentially random phase differences for 
different perturber configurations 0. In the following, 
spectral line broadening will be examined as an exam- 
ple of this behavior. For this case, two limits have been 
much studied. These are usually referred to as the Kubo 
fast(FML) and slow(SML) modulation limits. The tech- 
nique to be presented here, in fact, will be shown to be 
more widely applicable . 

The purpose of this letter is to propose a frequency 
separation for the treatment of the intermediate regime 
between the fast and slow modulation limits in the con- 
text of spectral line broadening in plasmas. Obtaining 
a theory that encompasses this intermediate regime is 
critically important in a number of physical problems, 
including Stark broadening, since the SML is attained ex- 
tremely slowly as the fluctuation frequency is decreased. 
With the frequency separation technique(FST) proposed 
here, practically all existing methods for the treatment 
of the intermediate regime overcome their problems. Of 
great practical importance for plasma spectroscopy is 
that the FST, in conjuction with the recently developed 
Frequency Fluctuation Model(FFM) HQ results in an 
ultrafast method that is valid for all parameter ranges 
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and all kinds of emitters. 

The FST idea comes from the following observation: 
The line profile is determined by V(t). V(t) is due to both 
electrons and ions, yet in line broadening calculations one 
almost always calculates their contributions separately. 
It must be emphasized that there is no mathematical 
justification for this step: The usual intuitive expla- 
nation is that for the vast majority of practical cases, 
either electrons or ions dominate the broadening, so that 
this separation should not introduce a large error ||. 
In fact, several tests of the electron-ion separation in a 
parameter range where electron or ion broadening was 
dominant have given excellent practical results. In the 
present letter we further give a general formal proof for 
this common belief by means of a Gedankencalculation. 

There remains the question of what happens when 
electron and ion broadening are comparable. There- 
fore to illustrate the separation hypothesis for compa- 
rable electron and ion broadening effects, a test cal- 
culation of the H a line emitted from a plasma with 
electron density of 4xl0 17 cm -3 and temperature of 1 
eV, has been performed. Fig.l presents the the dipole 
autocorrelation function, C(t), defined as the Fourier 
transform of the line profile, resulting from a calcu- 
lation of the emitting atom perturbed by a)electrons 
only(dashed-dotted), b)protons only (dotted) and c) elec- 
trons and protons (solid line). All calculations involve 
simulations performed with the non-interacting particle 
collision time-statistics method || using 1000 configura- 
tion averages, at which point convergence was achieved. 
The Gigosos et.al. method |l(| has been used for solving 
the Schrodinger equation and the code has been bench- 
marked against standard calculations P,[Tl||. The elec- 
tron and ion broadening contribute almost equally to the 
line shape as is evidenced by their individual C(t). The 
product of these individual correlation functions is rep- 
resented by the dashed line, and it can be seen that it 
is in good agreement with the correlation function cal- 
culated with both electrons and ions. This calculation 
indicates that the electron-ion separation is a valid as- 
sumption even in those rare cases where ions and elec- 
trons have comparable contributions, and where the sep- 
aration would seem to be least justified. 

Fig.l goes here 

In a unified treatment, that views electrons and ions on 
an equal footing, the main difference between the ionic 
and electronic fields is their fluctuation frequency. In 
this regard, the electron-ion separation is, in effect, a fre- 
quency separation of the total field into fast and slow(not 
necessarily quasistatic) components. 

There are two implications in the electron-ion separa- 
tion: First, to the extent that electrons may be treated 
in the impact approximation(IA), and to the extent that 
the IA is a solved problem, their part is solved and we 
only need at the end to convolve the ionic and electronic 
profiles. This is the main reason for neglecting electrons 
in our discussion and focusing on the ionic part from now 
on. In situations of practical interest, the ionic dynamical 
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behavior can range from the ion impact to the quasistatic 
regime and it is important to have methods capable of 
correctly treating ion dynamics. A particularly promis- 
ing modern method is the Frequency Fluctuation Model 
(FFM) plR]. This is an ultrafast stochastic method ca- 
pable of treating ion dynamical effects, whose usefulness 
for the analysis of todays highly complex experiments in- 
volving plasma spectroscopy cannot be overemphasized. 
It is in conjunction with the FFM that the frequency sep- 
aration proposed here should have its biggest impact in 
plasma spectroscopy. 

The second implication of the electron-ion separation 
is that if we were allowed to decompose the total electron- 
ion field into a fast electron part and a slow ionic part, 
can we do the same with the ion field alone? In other 
words, if Cif(t) and Ci S {t) denote the autocorrelation 
functions(AF) due to the interaction of the emitter with 
the fast and slow parts of the ionic field respectively(we 
specify later on what constitutes a fast or a slow compo- 
nent), and Ci(t) is the ionic AF, can one write 

Ci(t) w C if {t)C is (t) ? (2.1) 

We proceed to answer this question by means of the 
following Gedankencalculation: Consider an arbitrary 
stochastic interaction V(t), divided into fast and slow 
components and further suppose that for each particular 
realization (thereafter called configuration) of the inter- 
action, either the fast or slow component will be domi- 
nant. In such a case, suppose that in a large number N of 
configurations we have N s of them that are dominated by 
the slow component and N — N s dominated by the fast 
component. Let Yfc(i) denote the contribution to C(t) of 

the kth configuration, i.e. C(t) = J2k=i Y k(t)/ N - The 
weak component in configuration k will have ali(() ~ 1 . 
Consequently one obtains 

where C,(i) = Eji, Y{(t) and C 2 (t) = £Lv.+i Y' (i), 
with the superscipt s or f meaning that Y has been com- 
puted as a result of the action of the slow or fast compo- 
nent alone, respectively. Similarly, writing 

Cf(t)= Ns + E ^ Ykf{t) (2.3) 

and 

c M . »-«.+f~rt«) , (2 . 4) 

one may show that the absolute error, C(t) — C s (t)Cf(t) 
is [C s (t) — l][Cf(t) — 1]. This result justifies the sepa- 
ration approach in two cases: the case where the fast or 
slow component is dominant and also the case of short 
times, where C(t) is close to unity and the relative error 
is small, in agreement with the results of the numerical 
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simulations. The assumption in this derivation is that 
configurations where Yf and Y s are comparable are rare 
and that in the vast majority of configurations either the 
fast or slow component dominates. This is confirmed by 
detailed simulations and is also expected physically for 
any distribution function that is not too narrow. We fur- 
ther note that the argument given by Stamm, Smith and 
Talin that the slow component is essentially static 
during a fast collision and that one can neglect the en- 
tanglement of the fast and slow collisions applies equally 
to the fast-slow separation involving only the ion pertur- 
bation, as above. 

This Gedankencalculation is in fact quite general and 
applies both to the electron-ion separation and to the 
fast-slow separation of the ion component. In a prac- 
tical test (pi of the frequency separation on the ionic 
component alone, Eq.(l) was very well satisfied; how- 
ever, this test was done for a parameter range where the 
slow component was dominant. There remains the ques- 
tion of what happens when the fast and slow components 
are comparable, and this is a more important question 
in the context of the fast-slow ion field separation, be- 
cause there is no neutrality constraint relating the fast 
and slow components, as in the electron-ion separation. 
That is, it is clear that the FST is correct at short times 
and when either the fast or slow component is dominant; 
however, it is less clear when these components are com- 
parable. Although smoothness considerations imply that 
the error introduced by the FST should be small, this 
proposition must be confirmed by a calculation. Fig. 2 
illustrates the accuracy of the FST applied to the ion 
broadening in a parameter range where the fast and slow 
ion components are comparable, i.e. Lyman a transition 
from a plasma of singly ionized argon perturbers in a 
1.38 eV, 1.2x10 16 cto~ 3 plasma. Using the same meth- 
ods as for the calculation in Fig.l, we have calculated 
in Fig. 2 the autocorrelation function C(t) due to a) 
the fast ionic field component only (dotted line), b) the 
slow component only(dashed) and c) both fast and slow 
component (i.e. the net ionic) (solid line). Also shown 
is the product of a) and b) (dash-dotted line). All sim- 
ulation calculations involve 1000 configuration averages, 
at which point convergence was achieved. We note once 
again that Eq.(l) is very well satisfied. As a final remark, 
for this particular case electrons were neglected. Their 
broadening is almost an order of magnitude less but can, 
of course, be included by profile convolution. 

Fig. 2 goes here 

We have thus far said nothing on what constitutes 
a "fast" or "slow" component. The key point is that 
the fast component must be chosen so that it is treat- 
able by the IA. The IA is valid |li| for perturbers with 
HWHM, where v and p are the particle veloc- 
ities and impact parameters respectively and HWHM 
denotes the half width at half maximum. The way 
to "chop off" the fast component is then to determine 
a v/p = £1 > HWHM; perturbers involving higher 
frequency componcnts(i.c. higher v/p) belong then to 
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the fast component. Technical complications, associated 
with the precise meaning of the > in the above equation 
Q are deferred to a longer paper, as are practical ways 
to estimate the HWHM from the width of the impact 
component. 

At this point it should be stressed that we chose to 
separate out the impact fast component, not only be- 
cause the impact limit is a known limit, but also because 
slower, non-impact components are correlated and may 
not be frequency-separated: Only components varying 
rapidly on the memory loss time scale may be separated. 
An important implication of this fact is that at high den- 
sities a nonnegligible part of the electrons is no longer 
impact and becomes correlated with the ionic field. 

We now come to the implications of the frequency sep- 
aration for line broadening calculations: With the fre- 
quency separation the ionic calculation consists of com- 
puting the impact contribution(fast part) and the slow 
part contribution. It is clear that in this way the ion im- 
pact limit is recovered exactly and in as sophisticated 
a form as desired: For example, one may use the per- 
turbative IA |l3[ |, the exact semiclassical dipole results 
p5| (for hydrogenic emitters) or fully non-perturbative 
results including higher multipole effects |16fl . Further, 
this separation of the ionic field into fast(impact) and 
slow components improves the performance and/or reli- 
ability of every method capable of dealing with the in- 
termediate (ion dynamical) regime, which will now only 
have to correctly compute the contribution of the slow 
part: For simulations [1^], it is very expensive to have to 
integrate over sharp peaks. Such peaks are due to fast 
perturbers, which are correctly included in the impact 
part and do not appear in a simulation of the slow part. 
For the collective coordinate method ||, this separation 
is very important for the efficiency of the method, which 
would otherwise have to solve large linear systems. Fi- 
nally, stochastic methods have problems in attaining the 
ion impact limit. In the case of the -now obsolete- Model 
Microfield Method(MMM) @, which may only be valid 
for neutral hydrogen lines, the (perturbative dipole) im- 
pact limit is recovered with a delay, even for hydrogen. 
In the case of the FFM, the impact limit is not M a 
limiting case of the formalism, since it models the fluc- 
tuations of the interaction through a Stark component 
mixing model. Finally, experience with the BID |l9| ] 
method's approach to the impact limit is limited, but 
again the best that can be attained is the perturbative 
dipole impact limit. 

In the FFM (and also BID and MMM) all the time 
variation of the microfield is modelled via the covariancc 
of the electric field T EE (t) = (E(i)-E(O)). With the fre- 
quency separation, it is only the slow part of this covari- 
ance that will appear in the FFM(or BID); the fast part 
has been chosen so that it is treatable by the IA. 

As an example, the FST will be illustrated within the 
context of the FFM. The FFM will be used to describe 
only the slow component by employing a covariance with- 
out the rapidly fluctuating component. For example, to 
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describe the Tee(^) that represents only the slow compo- 
nent of the interaction between neutral emitters and in- 
dependent perturbers moving along straight trajectories 
||20| , the upper limit of the velocity integration can be re- 
placed by a finite value, f2/?, with p the impact parameter 
and f2 the frequency dividing the fast and slow compo- 
nents. This results in the expression for, \TEE(t)] a iow, 
the slow part of the ion field covariance: 

r sio 4i) = 2vm(-^) 2 / dv^[F(Q,b)-F(Q,-b)} 
47re J vt 

(2.5) 

where Q — t^t> ^ = vt/Xu and 



F(Q,b)= / dye- v e~-^ v2+b2+2b ^ y2 - Q2 (2.6) 



Q 

with \d the Debye length, n the density, f(v) the velocity 
distribution and Z the ion charge. 

It must be kept in mind that this result is given merely 
as an illustration, since more sophisticated models for 
^EE(t) may be used, either analytic [pi], or numerical. 
For example, if the covariance is determined from simu- 
lations, one should exclude the contribution to the total 
field of the particles giving rise to a fast component. 

Summing up, the frequency separation proposed in this 
work has solid theoretical as well as practical support and 
should prove very useful in improving existing methods 
dealing with the intermediate regime. The FST has been 
shown to work when either the fast or slow component 
is dominant and smoothness arguments, as well as prac- 
tical tests, indicate that it will indeed work for all cases 
of interest. The most important practical application is 
that with the frequency separation coupled to the FFM 
we have an ultrafast robust method, valid for all parame- 
ter ranges, a development of great importance for plasma 
spectroscopy. We note that the incorporation of the FST 
into the, now obsolete, MMM could furnish a unified pro- 
file for neutral emitters only, while the combination of the 
FFM and FST will provide a unified formulation for all 
emitters, neutral or charged for all plasma conditions. 
This now moves the focus of the outstanding problems 
to the computation of the strong interactions to the im- 
pact limit, which is a subject of much current interest 

Sim)- 
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III. FIGURE CAPTIONS 

FIG.l. Test of the independent electron-ion calcula- 
tion at a parameter range where electron and ion broad- 
ening are about equally important. 

FIG. 2. Test of the frequency separation on the ionic 
broadening. Particles with v/p > O = 1.44x 10 12 Hz 
have been included in the fast component and the rest 
in the slow component. The choice of f2 was made by 
equating the perturbative expression for the HWHM with 
one tenth of CI. 
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